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ABSTRACT 


This paper introduces the dual decomposition method for 
determining the distribution of an optimal objective function 
for a network problem. The objective function is to minimize 
the shortfall of demands to prioritized sinks for a four day 
period over a network that is subject to interdiction. The 
requirements of the model are that the upper and lower bounds 
on the capacities of the arcs and nodes of the network and the 
probabilities of interdiction are known. The dual 
decomposition method is an iterative approach to enumerating 
the possible instances of capacities in a capacitated network, 
based on the dual variables of the previous iteration. The 
purpose of the procedure is to determine the distribution of 
the shortfall of demands so that logistics planners can 
predict the performance of a supply distribution system over 


a short period of time. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. OBJECTIVE 

The objective of this paper is to provide a method for 
determining the approximate distribution of an optimal 
objective function of a linear program that minimizes the 
shortfall of supplies, over a period of several days, on a 
network that is interdictable. The supplies flow through arcs 
and nodes that have variable capacities to prioritized sinks. 
The procedure that is used to determine the distribution is 
the dual decomposition method. This procedure will be applied 
to an example network and the results will be analyzed. The 
result of the application of this method is that logistics and 
transportation Planners can predict the performance of the 
supply distribution system over a short period of time. 

The dual decomposition method is an iterative process for 
enumerating the possible flows through a network, based on the 
dual variables of the previous solution. This method is an 
attractive alternative in identifying the most likely patterns 
of performance of a network because complete enumeration is 
not required. The model is easy to use and understand and can 
be calibrated to meet the needs and capabilities of the 


planner. 


This method has obvious applications in the military 
logistics planning field. The logistics planners can use the 
method as a tool for determining the stockage levels of 
depots, the allocation of rear area protection forces to key 
segments of the transportation system, and to forecast periods 
of critical shortages. The predictive capability of the dual 
decomposition method allows the military commander to make 


more informed decisions on tactical matters. 


B. GENERAL 

The focus of this study is a military, land transportation 
network that has a limited capacity and is interdictable. The 
dual decomposition solution method can be applied to tactical 
and strategic movement routes that are multimodal, however for 
the purposes of this study the networks analyzed will be 
ground, tactical movement routes. The ground transportation 
system, whether it is highway, rail, or inland waterway, 
carries the vast majority of the required tonnages to the 
forward combat units. Therefore, it is the most important 
system to the military transportation planner. 

The transportation of supplies to military units involved 
in combat is a vital requirement of the United States Army's 
Air-Land Battle (ALB) fighting doctrine. The successful 


execution of this doctrine, as stated in the Army's Field 


Manual 100-5, will depend on four basic tenets: initiative, 
agility, depth, and synchronization. [Ref. 1] 

The intensity and speed of movements associated with the 
ALB doctrine will require the transportation movements planner 
and manager to be able to have fast, efficient tools to assist 
in the movement of supplies to support the combat units. The 
United States Army's Transportation Schools Field Circular 55- 
16, Movements Control Officer's Battle Support Guide, states 
that the planner must be able to maintain: 

- Uninterrupted support in spite of unit relocations; 

- Dispersion of supply activities; 

- The movement of cargo despite enemy interdiction 
throughout the theater and the depth of the battlefield. 
[Ref. 2] 

The requirements of the highly mobile, high-technology 
weapon systems will be enormous. The timely receipt of 
critical supplies will be vital to the successful execution 
of the ALB doctrine. The transportation community must be 
able to provide this uninterrupted support over networks that 
are subject to interdiction. 

Interdiction may be the result of sabotage by enemy 
SPETNATZ troops, enemy airstrikes, enemy missile strikes, and 
environmental factors, such as the weather. A high priority 
target of the enemy's assets will be the _  lines-of- 
communications (LOCs) and logistics facilities, such as depots 


and transportation centers. The securing of these objectives 


directly supports the Soviet concepts of surprise, mobility, 
deep penetration, and rapid exploitation. The Soviets are 
prepared to attack in considerable depth to destroy these 
targets and to harass the rear area operations (Ref. 3)-eethe 
transportation network may also be congested with the 
movements of friendly units and noncombatants which, in 
effect, interdict the network. 

This problem is not new to transportation planners. 
Considerable effort has been devoted to analyzing the problem 
of transporting supplies over a network that may be 
interdicted. The next section will review work that has been 
done in this area. 

Historical records indicate that military transportation 
planners have used extraordinary measures to ensure that 
critical supplies are delivered to units engaged in combat. 
For example, in World War II the Red Ball Express was 
developed to transport supplies to the combat units. The Red 
Ball Express was a term given to the military motor transport 
units that were centralized under a command with the mission 
of delivering critical supplies to the rapidly advancing 
combat units. The Red Ball Express was organized because the 
existing railways could not meet the demands generated by the 
armies' advance beyond the Rhine River. These units travelled 
on designated one-way roads so that there would be minimal 


interruption of supplies to the forward units. Despite its 


short life, the Red Ball Express was an effective method of 
delivering priority supplies. 

The Red Ball Express's accomplishments were impressive 
and widely acclaimed. During its eighty-one days of operation 
the transportation units carried an average of 5,088 tons per 
day. Despite the accomplishments of the Red Ball Express, the 
operation highlighted some serious deficiencies in the 
movement of cargo in a wartime theater. Centralized control 
of the movements was weak, the maintenance of the vehicles was 
poor, critical supplies, such as gasoline and tires, were 
consumed at an abnormally high pace, and it did not have 
sufficient cargo-handling capacity for long term operations. 
fmefi4) 

The Red Ball Express's average daily tonnage during its 
operation may have been sufficient in World War II, however 
today's modern weaponry and tactical systems require 
substantially more support. For example, a mechanized 
infantry division is expected to consume 20,000 tons of 
supplies during a day of offensive operations. 

It is necessary to capitalize on these lessons to prevent 
the lack of transportation support in future conflicts. In 
1946, Secretary of War Patterson said that in future wars the 
most efficient means of transport will be required and that 
",...changes in favor of speed and flexibility will make what 


we now have seem primitive." [Ref. 5] 


It is for this reason that military transportation 
planners must be able to analyze a transportation network in 
an effective and efficient manner so that the requirements 
are satisfied on time. By applying the dual decomposition 
method to a network the planner can determine the distribution 
of shortfall of supplies to the combat units. From this 
distribution the logistician and commander can analyze what 
effect this shortfall will have on combat operations, if any. 
This paper will illustrate one method that the planner can use 
to help ensure that the transportation system is receptive to 


the requirements of the modern battlefield. 


C. LITERATURE REVIEW 

The problem of analyzing the effects of interdicting 
networks has been studied for several years. There have been 
several different approaches to the analysis. One general 
approach has been to use linear programming and network flow 
theory and another approach has been to use stochastic methods 
and simulations to describe the performance of a network. When 
the arcs have variable capacities the different methodologies 
rely on the enumeration of all possible routes through a 
network. 

Ford and Fulkerson [Ref. 6] developed the maximum flow - 
minimum cut theorem. This theorem establishes that the 


maximum flow from a single source to a single sink is given 


by the value of the minimum cut that separates these nodes. 
This underlying theory forms a basis for developing more 
sophisticated network flow algorithms. 

Muster and McMasters [Ref. 7] developed an algorithm that 
determined which arcs in a network were the most susceptible 
to attack and how many of the available aircraft should be 
allocated to attacking a particular arc. The method used was 
a linear programming model that minimized the capacity of a 
cut set of a network subject to constraints on the number of 
available aircraft and the bounds on the capacity of the arcs 
of the network. The algorithm used a linear relationship 
between arc capacity and resource availability and was based 
on the maximum-flow minimum-cut algorithm developed by Ford 
and Fulkerson. Their method can only be used on planar 
networks that have a single source and single. sink. 
Additionally, the effect on the flow through the network is 
for one point in time rather than over a several day period. 

Preston and Howard [Ref. 8] developed a similar procedure 
to assign the optimum allocation of aircraft against a 
transportation network. Their method differed from Muster's 
and McMasters' in that they used an exponential relationship 
between arc capacity and resource availability in their 
algorithm. The purpose of the algorithm was to determine how 
to minimize the capacity of the network subject to aircraft 


availability. They determined the allocation of the aircraft 


against the network so that the incremental increase in number 
of aircraft to a mission is exceeded by the benefit resulting 
from disrupting the flow through the network. A limitation 
of this model is that the cost function is a measurement of 
the use of an aircraft versus the benefit of interdicting a 
network. The cost function must be varied in accordance with 
tactical considerations and is difficult to establish. This 
algorithm also has the same limitations as the Muster and 
McMasters paper. The algorithm is based on the minimization 
of flow for one day and the network has a single source and 
a single sink and was deterministic. 

Wollmer [Ref. 9] presented two algorithms for targeting 
strikes against a lines-of-communication network. His 
algorithms attempted to make the cost of achieving a 
circulation of flow between two points as large as possible 
over time. This was done by increasing the arc-cost functions 
and decreasing the arc capacities for a given period of time 
given the effect of targeting strikes. His first algorithm 
assumed the arc costs are linear functions of flow and his 
second one treated the arc costs as piecewise linear functions 
with one break point. The results of his algorithms enable 
the user to determine the effects of single and multiple 
strikes on the flow through a network for a point in time. 
Wollmer concluded that arcs that have a variable capacity 


based on the probability of interdiction can be replaced with 


arcs that have capacities that are equal to the expected value 
of the capacity. The conclusion that arcs with variable 
capacity can be replaced with arcs that have capacities equal 
to their expected value is incorrect. The expected value of 
an arc's capacity is a central value and the flow through a 
network is dependent upon the extreme values of the arc's 
capacities. Wollmer's algorithms can be applied to planar 
networks with one source and sink. 

Other authors used stochastic methods to analyze the 
effects of interdiction on networks. Four noteworthy papers 
are described below. 

Fishman [Ref. 10] described a Monte Carlo sampling plan 
for estimating the distribution of maximum flow in a directed 
network whose arcs have random capacities. His procedure 
determined an estimate of the performance of a network as time 
evolves given the joint distribution of the capacities of the 
arcs in the network. The network can represent a multistate 
system whose capacities are subject to random deterioration. 
While his paper does not deal with the network being 
interdicted directly, it does provide an estimate of the 
probability of the maximum possible flow through a network and 
the variance of this estimate. 

Bellovin [Ref. 11] analyzed a stochastic transportation 
network. He presented two related algorithms for analytically 


determining the value of the reliability of a network and the 


expected quantity of unsatisfied demand. The networks that 
were analyzed had random demands and the arcs and supply nodes 
were subject to failure with known distributions. Bellovin's 
algorithm computed the shortage distribution of a stochastic, 
small to moderately-sized, transportation network. His 
algorithm determined an expected quantity of the shortfall of 
supply rather than a distribution of the shortfall. 

Evans [Ref. 12] considered capacitated networks in which 
the branch capacity is an independent random variable with a 
known probability distribution. He investigated the 
computation of the probability distribution of the maximum 
flow in the network. His initial approach was to completely 
enumerate the capacity state vectors. The probability that 
the maximum flow has a specific value was found by summing 
the probabilities of all capacity vectors whose minimum cut 
set value equals the value of the maximum flow. Evans 
recognized that this approach produced computational 
difficulties for large networks. He then illustrated how a 
lattice structure can be used to reduce the computational 
requirements. The lattice is a set of cut sets of the 
network. Any two lattices that have comparable upper and lower 
bounds can be expressed as a single lattice. This approach 
may reduce some computations, but the number of cut sets 


grows exponentially as the network size increases. 
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Wollmer [Ref. i ls a) authored another paper on an 
interdiction model for sparsely traveled networks. He 
presented an algorithm for choosing locations to place forces 
in order to prevent the enemy from proceeding through the 
transportation network. His model calculated the probability 
of placing the force at particular arcs and nodes and the 
expected number of forces that should be placed on the 
network. His model used game theory as the approach of 
Placing protection forces against an infiltrator. This model 
did not deal with the degradation of arc capacities because 
of enemy attacks but it is useful in the assignment of 


protection forces on a transportation network. 


D. ORGANIZATION 

Chapter II presents the formulation of the four day 
logistics problem. The dual decomposition method will be 
detailed in Chapter III. Computational experience will be 
discussed in Chapter IV. The conclusions of the research and 
recommendations for further research will be addressed in the 


final chapter, Chapter V. 
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II. FORMULATION OF THE FOUR DAY LOGISTICS PROBLEM 


A. INTRODUCTION 

This chapter begins with an overview of network 
preliminaries. The formulation of the four day logistics 
problem for a combat setting and presentation of the model 
will be presented in subsequent sections. This chapter will 
provide the necessary background for the development of the 
dual decomposition solution method which is presented in 


Chapter III. 


B. DESCRIPTION OF NETWORKS 

A graph, G = (N,E), is a set of elements, N, and a set of 
unordered pairs of elements from N, E. The elements of the 
set N are referred to as vertices or nodes and the set E is 
the set of edges or arcs. The nodes are numbered by integers 
and are normally denoted by i and j. Nodes of a graph are 
numbered from 1 to n. Each edge is represented by (i,j), 
where i is the tail node (where the edge originates from) and 
j is the head node. 

A network is a graph with additional properties, 
parameters, or values associated with the graph's vertices 


and edges. 
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In a network representation of a transportation system, 
the nodes can represent supply depots, ports, transshipment 
points, critical road junctions, or a variety of other 
distinguishable features. The edges on the network represent 
the roads, rail lines, pipelines, or other avenues that 
supplies can travel over. Normally edges have labels that 
identify the edge's length, capacity, vulnerability, or any 
other characteristic or combination of characteristics. Nodes 
may also have properties and be labelled. 

The movement of supplies over a transportation network is 
analogous to a liquid moving through a pipe. The network, like 
a pipe, has a finite capacity. This capacity can be 
represented by such measures as vehicles per day, tons per 
Mile, or any other relationship which describes the limits of 
the arc's ability to carry goods. A network which has 
Capacitated arcs is called a capacitated network. The capacity 
of an arc may depend on the width of the road, the road's 
surface, weather, or the configuration and capacity of the 
vehicles that travel over the arc. 

Nodes may also have capacities. For example, a supply 
depot has a finite storage capacity. Again, these capacities 
can be described in much the same manner as an arc's 
capacities. 

Other characteristics of a network include sources and 


sinks. A source node is a node that originates flow and a 
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sink node is any node where flow terminates. In the 
transportation network, a source node may be a supply depot 
or port where cargo enters the network. A sink node may be 
a unit who consumes the supplies or an intermediate supply 
depot that stores supplies. Some algorithms only consider the 
network to have a single source and sink node. If the network 
being modelled has multiple sources or sinks and the objective 
function is the sum of flows, a super source and super sink 
can be constructed. A super source and sink can be 
constructed by adding dummy edges from the sink and source 
nodes to the respective super nodes. These dummy edges have 
infinite capacity. 

The underlying graph of the transportation network is 
assumed to be a sparse, planar, directed graph. A sparse 
graph is a graph where the cardinality of the edges, |E|, is 
proportional to the cardinality of the vertices, |N|. A graph 
is planar if it can be drawn in the plane with no arcs 
crossing. A graph is directed if there is only one 
orientation of the arc. That is, commodity flows can only 
move in one direction along the arc. 

The dual decomposition method does not require the graph 
be sparse or planar, however these attributes of the graph 
tell us about the performance of the method. 

The performance of a network can be measured and optimized 


by the application of linear programming methods. Some 
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reasonable objective functions of a linear program are to 
Maximize flow through a network, to minimize the cost 
associated with the movement of supplies through a network, 
or to minimize the time required to move supplies through a 
network. The objective function that is considered in this 
paper is to minimize the summed prioritized shortfall of 


demand at the sinks. 


C. MODEL DEFINITIONS 

The following section will define the terms associated 
with the model defined in the next section. 

The supplies that travel through an arc (i,j) during a 


time period, k, are indicated by f For the purposes of 


i-j-k* 
this paper, the supplies are the requirements of the combat 
units. These units require subsistence, petroleum and 
lubricants, barrier material, ammunition, repair parts, and 
major end items such as trucks and tanks. 

If a demand is not satisfied at a sink a shortfall will 
occur. The shortfall of supplies at a sink j during a time 
period, k, is represented by sf;,. The supplies that do get 
to a sink plus the shortfall equal the demand. 

The maximum capacity of an are that has not been 
interdicted during a specific time period, k, is represented 


by Teas Each arc is vulnerable to interdiction by the 


enemy. The vulnerability of the arc can be expressed as a 


aS 


probability that the arc will be interdicted and therefore, 
have a reduced capacity, ay es The probability that arc 
(i,j) is interdicted on day k is represented bY Pi.j.K- The 
probability an are is not interdicted is 1 - Dijk - ine 
capacity of an arc depends on the width of the road, the 
road's surface, weather, the capacity of vehicles using the 
arc, and whether the arc has been interdicted. 

The demand at a sink node, j, for a time period, k,is 
indicated by D;,- This demand, measured in short tons and 
gallons, can be determined by summing the requirements for 
all commodities for that particular node or by standard 
planning factors. Standard planning factors are guidelines 
that provide the planner some rules of thumb on the expected 
consumption rates of different commodities. They are based on 
historical usage rates of the different supply commodities. 
They can be used to estimate the requirements of a particular 
type of unit during the planning phase or when the 
requirements are not explicitly known. These factors are 
available in U.S. Army Field Manual 101-10-1, Staff Officers' 
Field Manual, Organizational, Technical, and Logistic Data, 
and Reference Book 101-999, Staff Officers' Handbook. It will 
be assumed that the supplies entering the network are not 
constrained by availability. 

The inventory of supplies that is held at a depot node, 


j, between time periods k-1 and k is denoted by Ty. The 
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capacity of a node, j, for a time period, k , is denoted by 
CiiKe The capacity of a node can be restricted by the size of 
the facility, availability of material handling equipment, 
and a variety of other factors. The sum of the flows into a 
node minus the sum of flows out must be less than the capacity 
of the node during time period k. If the node is at its 
maximum capacity and is not being interdicted the capacity is 
denoted by Cai The notation for a depot that is being 
interdicted is indicated by c’,,. 

A weighting factor for a particular sink node and time 
period is represented by w,;,.- The weighting factor allows 
the planner to prioritize the shipment of cargo to a 
particular sink in accordance with the priority of support 
determined by the commander. A large weighting factor implies 
the node has priority and shipments should arrive in a timely 
manner. Each sink has a different priority, weighting factor, 
that corresponds to the priority established by the commander 
of the forces receiving the supplies. 

The set of nodes that represent the sinks is denoted by 
T. The set Dis the set of all depot nodes. The problem will 
be solved for K days. 

The objective function of the model is to minimize the 
summed, weighted shortfall of the demand for each sink and 


time period. 
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The variables used in the model are summarized below in 


Table 1. 


TABLE 1 SUMMARY OF LINEAR PROGRAM VARIABLES 


fij. 7 flow from arc i to arc j on day k 
sii. - shortfall of supplies at node j during day k 
Vereen - the maximum capacity of an uninterdicted arc (i,j) 


during day k 


a oe the capacity of an interdicted arc (i,j) during 
day k 

Pi.j.x 7 the probability that an are (i,j) is interdicted 
on day k 

D; x - the demand at node j during day k 

1 55 - the inventory at node j between day k-1 and k 

c’,, - the uninterdicted capacity of node j during day k 

(oe Beh - the interdicted capacity of node j during day k 

Wk - priority of sink j on day k 


D. FORMULATION AS A LINEAR PROGRAM 

The problem for the logistician and the movements planner 
is to determine the distribution of the shortfall of supplies 
so that some insight can be gained on the performance of the 
Supply distribution system over a period of time. The 
distribution of the optimal shortfall involves minimizing the 
shortfall of supplies to the prioritized sinks for each day. 
The minimization is subject to the balance of flow in and out 


of the depot nodes and the conservation of flow at non-depot 
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nodes. Additionally, constraints on the capacities of the 
arcs and depot nodes exist. 


The problem can be written mathematically as: 


Minimize >" Wie (SEj,) 1b eal 
j€sinks k-1,2....., K 
Subject to 
het st, = Dj. for all je€T de 
ié€Nn k = V2 fone pK 


- poe fees Pog «=O COO for all j € D 153 
iéN k = Lis Sai ae 4 
ee eet lias ae (Wao le Fee Wee! 
i€N ié€N 
for all j € D 
koe" 2534.5 +, aK 
i ee ic 5 Dey a cea 1.5 
i€Nn i€N 
for all 4 N - {D U T} 
for ke= 1 Ze. sak 
O < fy k < Ui jk for all aA) € N 1.6 
k ow, 2....,K 
0 Sle Ca £or ail 3 € D Ve) 


Jr ai Ue Ae ene 44 
Equation 1.1 is the objective function. The objective is 
to minimize the shortfall of supplies at sinks according to 


a certain priority or weighting factor for four time periods. 
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Equation 1.2 is the constraint determining shortfall from 
flows into the sink node. It states that the flow of supplies 
into a sink plus the shortfall equals the demand. 

Equation 1.3 states that the flow of supplies out of a 
depot cannot be greater that the inventory of the depot. This 
equation, and Equation 1.4, imply that supplies cannot move 
directly through a depot node in one time period. These 
equations ensure a proper time sequencing of the flow of 
supplies through the network. A requirement of this model is 


that depot nodes must not be greater than one time period away 


from each other. The model is not limited by this 
requirement, however. The time period can be defined in 
whatever time increments that are suitable to the user. For 


this paper the time periods will be in terms of days. 

Equation 1.4 states that the flow of supplies into a depot 
plus the previous day's inventory equals the flow out of a 
depot plus the ending inventory for that day. 

Equation 1.5 is a constraint on the non-depots. It states 
the flow of supplies into a node during a time period must 
equal the flow of supplies out of a node during the same time 
period. This ensures that supplies are not stored at non- 
depot nodes. 

Equations 1.6 and 1.7 are restrictions on the arcs and 
depots, respectively. The flow of supplies through an arc 


and node must be less than its maximum capacity but larger 
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than zero. The constraint on the lower bound of the 
capacities ensures that arcs are one directional and that 


depots can only ship supplies that are on hand. 


E. NETWORK MODELLING TECHNIQUES 

The linear programming problem described in the last 
section can be transformed to a network flow problem. This 
transformation enables the problem to be solved using more 
efficient and faster network solvers such as GNET [Ref. 14]. 
A network flow problem also produces solutions that are 
integer. Another benefit of the transformation to a network 
structure is that the physical characteristics of the problem 
are better represented by a network flow problem and they are 
more readily accepted by nonanalysts [Ref. 14: p. 2]. 

Some modelling techniques are required to transform the 
problem from a linear programming model to a network flow 
model. Artificial arcs are used to represent the flow of the 
inventory of supplies stored at depot nodes from one day to 
another. Artificial arcs are also required to maintain the 
conservation of supply and demand. In a network model the 
tonnage demanded must equal the tonnage supplied at the supply 
node. In the linear programming model described in the last 
section, a shortfall can exist at a sink. The objective 
function value is the amount of the shortfall times a priority 


factor. Since the objective of the network model is the same, 
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an artificial arc is used to capture the shortfall. In a 
minimum cost network flow problem the shortfall can be 
determined by having costs associated only with the arcs that 
carry the shortfall. 

Another technique that is required for the transformation 
to a network flow model is the addition of another arc and 
node for every depot node. The depot node in the linear 
programming model has capacity constraints. The network model 
also requires constraints on the capacity of the depot 
inventory but only arcs in the network formulation can be 
constrained. It is necessary to construct an arc that 
represents the capacity of the depot node and another node to 
link the arc to the remainder of the network. This arc does 
not have the same physical interpretation as the other arcs 
in the network, but it does allow the capacity of the depot 
node to be constrained. 

The next chapter discusses the dual decomposition method 


and its application to a network that can be interdicted. 
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III. THE DUAL DECOMPOSITION METHOD 


A. INTRODUCTION 

This chapter gives the conditions under which the dual 
decomposition method can be applied to a network flow problem 
and the steps of the dual decomposition method. Then the dual 
decomposition method is applied to a simple example to 


illustrate how the methodology works. 


B. APPLICATION OF METHODOLOGY 

The dual decomposition method can be applied to a network 
where the arcs have known higher and lower capacities. The 
higher capacity is the capacity of the arc when it is not 
interdicted and the lower capacity is the capacity of the arc 
when it has been interdicted. Generally, the capacities are 
expressed in terms of rates of flow, such as tons per day or 
gallons per hour. The values of these capacities can be 
determined from historical data, subject matter experts, and 
some capacity values can be found in the DAMSEL network data 
base [Ref. 15]. This data base has unrestricted capacity 
values for the majority of the transportation networks in 
Europe. A lower capacity of an arc may be zero but it has 
been determined from studies of network capacities in wartime 


that an arc in the network normally has some positive value 
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because the damage done by interdiction can be repaired in a 
short amount of time [Ref. 16]. 

The probability of interdiction must also be known. 
Again, the probability of interdiction can be determined in 
the same manner and from the same sources as the arc 
capacities. The probability of interdiction normally depends 
on how the enemy perceives the tactical or strategic value of 
the arc. 

The dual decomposition method determines the distribution 
of the optimal objective function of the minimum cost network 
flow problem without having to enumerate all possible 
combinations of arc capacities. This is especially critical 
given that each arc has two possible arc capacities and there 
may be hundreds of arcs in a network. If there are n arcs in 
a network, the determination of the true distribution of the 
optimal objective function by full enumeration will require 
solving 2” minimum cost network flow problems. It will be 
illustrated that the dual decomposition method can determine 
an approximate optimal objective function with substantially 
fewer problems being solved. 

In deterministic optimization, the expected value of the 
Capacity of each arc is used. A single minimum cost flow 
problem is then solved. The dual decomposition method also 
precludes having to use an expected value for the arc's 


Capacity. The use of expected values for arc capacities does 
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not provide an expected value of a minimum cost network flow 
problem because the minimum cost flow depends on extreme 
values of the arc capacities rather than central values. This 
will be illustrated during the example problem later in the 


chapter. 
C. DUALITY AND THE MULTIDAY LOGISTICS PROBLEM 


The multiday logistics problem presented in Chapter II may 
be rewritten as: 
Minimize c’f = 2Z crue 
Subject to Af =b 
Tf <s u 
£ > Oo 
The vector of flows is represented by f, c is the cost vector, 
and A is a K(!N! + |DI) by KIEI matrix encapsulating the 
Sonstraint equations 1.1 through 1.5. I is a K(ID!| + IEl) 
dimensional identity matrix. The vector of arc capacities is 
denoted by u. 

The dual to the above linear program is given by: 
Maximize b’v + u'w = Z Bee 
Subject to vA+ wI > ¢c 

w> 0 
The vectors v and w are the dual variables of the problem. 
A basic result in dualtiy theory is that, at optimality, the 


objective function values of these two problems are equal. 
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Let it E, where E is the set of edges, then dZ/du, = wj, 
as evidenced by the objective function of 3.2. Hence, for a 
given set of capacities the arc with the largest optimal dual 
value can be seen to be the arc which is the most crucial in 
determining the value of the objective function, Z. 

The network solver GNET was used to solve instances of the 
minimum cost network flow problem represented by the linear 
program given at 3.1. The GNET code was modified in order to 
determine the vector w post optimally. This modification does 


not significantly impact the performance of the GNET code. 


D. THE DUAL DECOMPOSITION METHODOLOGY 

As has been mentioned, we are interested in developing an 
iterative approach to the decomposition of the stochastic 
minimum cost network flow problem. The methodology involves 
starting with all the arcs "free", that is, each arc may be 
at its upper or lower capacity, the capacity being still 
random. The set of free arcs will be referred to as the set 
12 5 

The starting problem is called "problem 0." At the 
outset, Fy = E. The subscript indicates the problem number. 
Two minimum cost network flow problems are solved, subproblem 
A has all are capacities set to their lower bounds, while 
subproblem B has all arc capacities set at their upper bounds. 
After solving these two problems using the modified version 


of GNET, the dual vectors w, and wz, aS well as the objective 
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function values Z, and Z,, are determined. Recall from the 
discussion of the dual of the multiday logistics problem that 
for each arc i, wa, is the change in the objective function 
value for one unit of change in the capacity of arc i given 
that all other arcs stay at their lower capacities. Thus the 
arc i* which maximizes Wai; is the arc of greatest influence in 
Subproblem A, while the arc 5. which maximizes Wg,j iS the are 
of greatest influence in subproblem B. Among A and 3) arc 
k" is chosen which maximizes {w,, : i¢€ E} U {wg,; : J€ E}. Are 
k" will be called the "pivot arc" for problem 0. tie ae 
condition on the capacity of are k", we will gain significant 
insight into the value of the objective funtion. 

Let problem 01 have the capacity of arc k set to its 
upper capacity, while problem OO has arc k*" set to its lower 
capacity. The free arcs for both problems are given by Fj) = 
Fo) = E - { k }. We assign probabilities of occurrance to 
each of these new problems. Pp), is the probability that arc 
k” is set to its upper capacity. In the original problem 
formulation, this corresponds to arc k* being unharrassed. Pg 
is given by the probability that arc k is interdicted, 
implying its lower capacity prevails. 

The next step in the process is finding values of Z, and 
Z, for problems 01 and OO. For problem OO, both subproblems 
have the capacity of arc ko (all the arcs in E - F) set to its 


upper capacity. In problem 01 subproblem B, the free arcs are 
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set to their upper capacity to find the new Z, value. In 
problem O01 subproblem A, the free arcs are set to their lower 
capacity. Note that, the new Z, value equals the Z, value of 
problem O because the arc capacities for all arcs are the 
same. 

In problem 00 subproblem B, the free arcs are set to their 
upper capacity to find the new Z, value, for subproblem A of 
problem 00, the Z, value is found by setting the free arcs to 
their lower capacity. Every arc in E - F, that is arc k* in 
this case, is set to its lower capacity. The new Z, value 
equals the Z, value of problem O because, again, the arc 
capacities are identical. Thus in the first pair of problems, 
we have only needed to solve two new subproblems. 

In problem O1, the vectors of dual variables for both 
subproblems are analyzed to find the new arc k* which 
maximizes {w,; : ie E} U {wg,; : Jj € E}. This gives the pivot 
arc for problem Ol. 

In problem 00, the same procedure is used. The vector of 
dual variables is analyzed to find the arc k* which maximizes 
(Wa; : i¢€ E} U {wy ; : j€ E} for problem 00. This yields the 
pivot arc for problem OO. 

We continue by generating problems 001 and O00 from 
problem OO and problems O10 and O11 from problem 01. 

The next step in the dual decomposition process is to 


check the termination criteria for both problems OO and Ol. 
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A problem can be terminated if it can be determined that 
further decomposition of the problem will not produce 
significant improvements in the approximation of the 
distribution of the optimal objective function. The 
termination criteria measures the absolute difference between 
the problem's Z, and Z, values multiplied by the probability 
that the nonfree arcs realized their prescribed values for 
that problem. The equation for the termination check is 
provided at Equation 3.3. 

IZ, - Z| * p < threshold value 3.3 
Where p = Po, for problem 00 and p = Py, for problem O01. 

The threshold value is determined by the user of the 
methodology. The threshold value must be small enough so that 
the approximate distribution of the optimal objective function 
that results from the application of the methodology 
adequately represents the actual optimal objective function. 

The iterative process of constructing subproblems A and 
B yielding Z, and Z,, respectively, checking for termination, 
finding the pivot arc and probability, and creating the two 
daughter problems continues until all problems terminate. 

If the termination criteria is met, the values of Z,, Z3, 
and p are added to a SOLN set. The SOLN set is the set of 2Z;, 
Z,, and p values that will determine the approximate 


distribution of the optimal objective function. 
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The probability density function of the minimum cost flow 
for the network may be constructed by plotting the values of 
Z,, Z2, and p from the SOLN set. For each terminated problem 
we assume that the objective function value for a terminated 
problem is uniformly distributed between Z, and Z, with 
probability p. The result of the combination of the 
individual Z,, Z,, and p values is the approximate distribution 


of the optimal objective function. 


E. EXAMPLE PROBLEM 

This section contains an example problem that illustrates 
the dual decomposition method. For the purpose of 
illustrating the methodology, a different, simpler model than 
the four day logistics problem described in Chapter II will 
be used. This linear program's objective is to maximize the 
flow through the network. The only constraints are the 
balance of flow constraints through the nodes and the capacity 
constraints of the = arcs. The model is described 
mathematically below. 


maximize f 


such that 
f ase 3k Sh 
= a) aa 0 if i*«lor 7 
je Nn m€N 
-f guse SL = We 
iy Sf Ci; 
> 0 


Xi. j 
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The flow through the network is represented by f, the flow 
through an are (i,j) is represented by x;,,, and the upper 
capacity of are (i,j) is represented by Ci 4° 

The example is a simple network with four arcs and four 


nodes. The example network is represented at Figure 1. 


Source 


(8,6,0.3) (5,4,0.7) 


ARC B ARC D 


(high cap, low cap, prob of interdiction) 





Figure 1. Example Network 


The threshold value of the termination criteria will 
arbitrarily be set to 0.3 for this example. 
The first step in the method is to solve the problem when 


the arcs are at their lower capacities and find the objective 
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value, Z,. The problem is then solved at the upper capacities 
to find Z,. When all arcs are at their lower capacities, the 
flow through the network, Z,, is 7 units because the flow is 
restricted to the maximum capacity of Arcs A and D. The value 
of Z, is 10 because Arcs A and D restrict the flow to 5 units 
in each branch of the network. The probability, p, is 1 since 
there is no pivot for the initial problem. 

The termination criteria is checked and since the value 
is above 0.3 the decomposition begins. 

The next step is to examine the dual variables of the two 
solutions and determine which arc has the most negative 
reduced cost. The arc with greatest impact on Z, and Z, values 
is Arc A. 

The decomposition process begins with the problem being 
pivoted on Arc A. Arc A is fixed at its lower capacity and 
the problem is solved with the remaining free arcs at their 
lower capacities and then again at their upper capacities. 
The value of the objective function for the arcs at the lower 
capacity is the same as Z, of the previous problem. The value 
of Z, for this branch of the problem is 8. 

Arc A is then fixed at its upper capacity and the problem 
is solved with the remaining free arcs at their upper and 
lower capacities. The value of Z, in this case is the same as 
Z, of the previous iteration. The value of Z, for this branch 


of the decomposition is 8 (see Figure 2). 
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The termination criterion for both problem sets is checked 
next. The value of the termination criterion for the problem 
set with Arc A set its upper capacity is f8- => 40! * 0.6 = 
1.2. The value of the termination criterion for the problem 
set with Arc A set at its lower capacity is 0.4. Since both 
problem set's solutions do not meet the termination criterion, 


the decomposition continues on both branches of the problem. 


72 = 10 
jz1 - 22) *P = 


Dies 
Tea No 
ieeZ2 |e |= 3 Z1 
22 
Tio7.2\ *“ira=-0.4 


(3,0.4) 





Figure 2. First Branch Diagram of Example Problem 
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The dual variables for the problem pair that had Arc A 
set to its upper limit indicated that the decomposition should 
continue with Arc C being fixed. The decomposition generates 
two additional problem pairs. The first problem pair is with 
the upper capacity of Arc C fixed and the values of Z, and Z, 
determined for the new branch. As in the first decomposition, 
the value of Z, for this branch of the solution equals the 
value of Z, from the previous solution, which was 10. The new 
value of Z, for this problem set is 9. 

The second problem pair has the lower capacity of Arc C 
fixed and the new value of Z, is determined. The value of Z, 
for this branch equals Z, of the previous solution. The new 
Z, value is 8 and the value of Z, from the previous solution 
was 8 also. 

The termination criteria values for both problem pairs 
are then calculated. The first problem pair had a termination 
criterion value of 0.12 and the second problem pair had a 
value of 0.0. Since both of these values are less than the 
threshold value, the decomposition on this branch of the 
problem is stopped. The values of Z,, Z2, and p for both 
problems in this pair are added to the SOLN set. Figure 3 is 
a branch diagram of the problem at this point in the 


decomposition. 
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Figure 3. Second Branch Diagram of Example Problem 


The analysis of the dual variables for the problem pair 
that had Arc A set to its lower capacity determined that Arc 
D should be pivoted on. The decomposition of the network 
continues at this branch with Arc D being set to its upper 
and lower capacity and two new problem pairs being generated. 
These two new problem pairs are solved and the termination 
criteria checked. If the termination criteria are not met the 
decomposition process continues. 

Since both branches of the problem met the termination 
criteria the decomposition process ends. A branch diagram of 
the final structure of the problem, with all remaining values 


provided, is at Figure 4. 


22 





Figure 4. Final Branch Diagram of Example Problem 


The distribution of the optimal objective function is then 
determined. The values in the SOLN set are plotted on an X-Y 
axis. The X axis is the values of the network flows, Z, and 
Z,, and the Y axis is the conditional probability. The Z, and 
Z, values determine the upper and lower limits on the network 
flow for that particular case. The probability density 
function of the network flow for the example problem is 


provided at Figure 5. 
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Approximate Optimal Objective Function 
for Example Problem 


Probability 








Figure 5. Approximate Optimal Objective Function 
for Example Problem 
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F. COMPARISON OF RESULTS OF DUAL DECOMPOSITION METHOD AND 
OTHER METHODOLOGIES. 


Determining the exact distribution of the flow for the 
example network would require the complete enumeration of all 
possible flow combinations, solving 24 = 16 minimum cost 
network flow problems. The dual decomposition method only 
required solving 6 minimum cost network flow problems. 

The solution of the network flow problem by using expected 
values for arc capacities only requires solving one network 
flow problem, however the result is not an accurate 
representation of the actual flow through the network. 

Table 2 provides a summary of the comparison of the 
different results. 


TABLE 2. COMPARISON OF RESULTS OF EXAMPLE PROBLEM 
BY DIFFERENT METHODOLOGIES 





TABLE 2 


COMPARISON OF RESULTS OF EXAMPLE PROBLEM BY 
DIFFERENT METHODOLOGIES 


Type of Methodology Flow Problems Required 


Full enumeration 8.02 
(actual) 


Dual Decomposition 7.9 


Expected Value 8.7 





The result of the problem by full enumeration is the 


distribution of the actual flow through the network, however 


38 


it requires solving 10 additional problems. The dual 
decomposition method provides a suitable approximation to the 
optimal objective function and only requires solving 6 
problems. The solution by the expected value procedure is an 
overestimation by 8.7% of the true flow and does not provide 
an accurate representation of the flow through the network. 
The variance of the full enumeration solution is 0.807. 
Therefore, the solution determined by the expected value 
method is not within one standard deviation of the actual 
solution. The solution determined by the dual decomposition 


method is well within one standard deviation. 
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IV. COMPUTATIONAL EXPERIENCE 


A. INTRODUCTION 

The first section of this chapter presents the results of 
a four day logistics problem when the dual decomposition 
method is applied to a network representation of the problem. 
The second section examines the effects on the approximate 
distribution of the optimal objective function when the 
threshold value of the termination criterion is varied. The 
third section presents the structure of a computer program 
that uses the dual decomposition method to determine the 
approximate optimal objective function for a minimum cost 
network flow problem. 
B. RESULTS OF A FOUR DAY LOGISTICS PROBLEM DETERMINED BY THE 

DUAL DECOMPOSITION METHOD 

A logistics scenario was developed to determine the 
approximate shortfall of supplies to two combat units. The 
linear program that represented the scenario was transformed 
into a minimum cost network flow model by using the techniques 
described in Chapter II. 

The network representation of the logistics problem 
consisted of 57 arcs and 48 nodes. The objective of the 
problem was to minimize the shortfall of demands to two combat 


units over a four day period. Each unit had a priority of 
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support that varied with each day. Each arc had a higher 
Capacity, a lower capacity, and a probability of interdiction. 
The network had 12 depots, 12 non-depots, 8 sinks, and one 
source. The network represented a fictitious highway network 
with capacities, probabilities, priorities, and demands 
established arbitrarily. The two combat units had a combined 
demand of 3125 units of prioritized supplies during the four 
day period. 

The dual decomposition method was used to determine the 
approximate optimal objective function of the minimum cost 
network flow. The initial threshold value for the termination 
criterion was 7.0. The determination of the approximate 
distribution of the shortfall of supplies required solving 64 
minimum cost network flow problems. After the termination 
criterion had been met for each pair of problems, the 
approximate distribution of the optimal objective function of 
the problem was constructed. The probability distribution 
function of the optimal objective function is in Figure 6. 

It was determined from the probability distribution 
function that the prioritized shortfall of supplies to the 
two combat units for the four day period would have extreme 
values of 1350 and 1730 units of prioritized supplies. The 
mode of the approximate distribution was at 1550 units of 


prioritized supply. This is the most probable value of the 
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Approximate Optimal Objective Function 
for a Minimum Cost Network Flow Problem 
(Threshold Value = 7.0) 


Probability 


1450 1550 1650 
Shortfall of prioritized supplies 








Figure 6. Approximate Optimal Objective Function for a Minimum 
Cost Network Flow Problem (Threshold Value = 7.0) 
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shortfall of prioritized demand. There are several "spikes" 
in the distribution, or local modes, between 1600 and 1680 
units. The logistics planner can approximate that the 
shortfall of demand will be between 1550 and 1650 units in 
the shortfall of supplies 80% of the time. 

Since the shortfall of supplies would be approximately 
50% of the requirements, the logistics planner could recommend 
to the combat units’ commander alternative methods of 
distributing the supplies and possibly how the shortfall could 
be avoided by increasing the stockage levels of supplies at 
the forward depots. The logistics planner could also 
recommend that certain segments, arcs, in the transportation 
network receive additional protection from interdiction 
thereby decreasing the probability of interdiction. 

If the problem described above had been solved by complete 
enumeration of the arc capacities, the actual shortfall of 
supplies could be determined. However, this would require 
Bo) (eve approximately 1.4 x 1077, minimum cost network flow 
calculations. The number of solutions required by the dual 
decomposition method is 2°) = 4.44 x 107° of the total number 
of solutions required by full enumeration. 

The logistics planner normally only requires -= an 
approximate value of the shortfall of supplies. The dual 


decomposition method can accomplish this approximation by 
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solving a very small fraction of the calculations that are 


required by the full enumeration procedure. 


C. EFFECTS OF VARYING THE VALUE OF THE THRESHOLD VALUE 

If the threshold value of the termination criterion is 
adjusted to 15.0, only 32 minimum cost network flow problems 
must be calculated. This is a 50% reduction in the number of 
calculations required when the threshold value is 7.0. 

All branches of the final decomposition tree were the same 
length. This will not normally be the case because the branch 
with a low probability will satisfy the threshold criterion 
earlier than the branch with a higher probability. In the 
four day logistics problem the probability of interdiction for 
each arc was normally 0.5 or 0.4. The probabilities for each 
branch in the tree were similar so the threshold criterion was 
satisfied at the same stage in the partial enumeration. 

The approximate distribution of the optimal objective 
function for the minimum cost network flow with a threshold 
value of 15.0 is at Figure 7. 

This distribution also has a mode at 1550 units of 
prioritized supplies but has a relative mode at 1500 units. 
The approximate distribution with the threshold value set to 
15.0 is smoother than when it is set to 7.0. The logistics 
Planner can approximate that 80% of the shortfall of demand 


will occur between 1490 and 1650 units of prioritized 
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Approximate Optimal Objective Function 
for a Minimum Cost Network Flow Problem 
(Threshold Value = 15.0) 


Probability 


| | r+ 
1460 1860 


Shortfall of prioritized supplies 





Figure 7. Approximate Optimal Objective Function for a Minimum 
Cost Network Flow Problem (Threshold Value = 15.0) 
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supplies. This 80% confidence interval for the approximate 
distribution is wider than when the threshold value is 7.0 
but it still can be used to determine an approximate shortfall 
of supplies to the combat units. 

The following table illustrates the results of the four 
day logistics problem when it is solved using the dual 
decomposition method and when expected values are used for 
arc capacities. 


TABLE 3. COMPARISONS OF RESULTS OF FOUR DAY LOGISTICS 
PROBLEM SOLVED BY DIFFERENT METHODOLOGIES 


TABLE 3 


COMPARISON OF RESULTS OF FOUR DAY LOGISTICS PROBLEM 
SOLVED BY DIFFERENT METHODOLOGIES 


Dual Decompositton Method Are Expected Value 


Threshold Value 
7.0 15.0 


Mode 1660 1650 


Expected 1666 1570 
Value 


Std 
Deviation 





The results illustrate that the solution of the method 
that uses expected values for arc capacities overestimates 
the capabilities of the network. The solution for this method 


has a lower value for the shortfall of prioritized supplies 
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than the solutions that used the dual decomposition method. 
This lower value for the shortfall of supplies gives an 
inaccurate and more optimistic prediction on the network's 


capabilities. 


D. DESCRIPTION OF THE COMPUTER PROGRAM 
The steps in a computer program for solving the multiple 
day logistics problem are provided below. The algorithm is 
written in pseudo-code so that it can be implemented in any 
programming language. 
Read in the network characteristics 


Construct the multiple day network flow representation of 
the problem 


Construct the head, tail, capacity, cost, and probability 
arrays 


Solve for Z, 
Solve for Z, 
IF( {Z, - Z,| * P < threshold value) then 

Record Z,, Z,, and p values 

STOP 
ELSE 

100 

Choose arc i to be pivoted on 

Put a new problem on the queue with arc i set to 

its upper capacity 
Put a new problem on the queue with arc i set to 


its lower capacity 
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Take a problem off the queue 
Set all free arcs to upper capacity 
Call GNET 
Record primal and dual variables 
Set all free arcs to lower capacity 
Call GNET | 
Record primal and dual variables 
IF( |Z, - Z,| * P < threshold value ) then 
Record Z,, Zz, and p values 
IF( the queue is empty ) then 
STOP 
ELSE 


GO TO 100 


The queue that contains the problems to be solved is a 
first-in-first-out queue. The array that represents the 
problem to be put on the queue has values of Z, or Zz, and p 
from the parent problem pair. If the problem being put on 
the queue is going to determine a new minimum cost network 
flow for the arc set to its lower capacity, the Z, value is 
recorded. This prevents having to solve for the Z, value, 
for that problem pair, since it will be the same as the last 
problem pair. The opposite holds if the arc is set to its 


upper capacity. The Z, value is recorded since it is the same 
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as the previous problem pair. The conditional probability for 
the previous problem pair is also recorded since it will be 
used to determine the new conditional probability. The array 
that is stored on the queue also has a dimension that contains 
a cell for every arc. If the arc is not free the 
corresponding cell for that arc is marked so that when the 
problem is taken off the queue, the arcs that are not free can 
be set to their upper or lower capacity. The arcs will be set 
to their upper or lower capacity depending on how they are 
marked in the cell. 

A FORTRAN program that uses the dual decomposition method 
to determine an approximate distribution of the optimal 
objective function for a minimum cost network flow problem is 
provided at Appendix A. This program uses the network solver 
GNET to solve the minimum cost network flow problem. The 
program also constructs a multiday network representation of 
the problem given the road network and depots. 

The next chapter provides conclusions on the use of this 
methodology on network interdiction models and then provides 
recommendations on areas of further research that should be 


explored. 
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V. SUMMARY AND FUTURE RESEARCH AREAS 


A. SUMMARY 

The dual decomposition method is an attractive alternative 
to solving minimum cost network flow problems when an 
approximation of the distribution of the optimal objective 
function can be used. This method can determine a close 
approximation of the distribution of the optimal objective 
function by solving a fraction of the number of problems 
required by full enumeration. 

A close approximation of the distribution of the optimal 
objective function is adequate for planning purposes in most 
cases. The dual decomposition method allows this 
approximation to be determined by solving fewer problems then 
other methods, so it is more responsive to the time 
constraints of the logistics planner. 

The dual decomposition method is also preferred over using 
expected values for arc capacity values. The expected value 
approach overestimates the optimal objective function and does 
not provide an accurate prediction on the capability of a 
transportation network. This overestimation on capability 
could produce unacceptable consequences if predictions on the 
logistics support available to units engaged in combat are 


based on the results of this method. 
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The dual decomposition method can also be used to make 
policy decisions. If the results of the analysis of a 
transportation network by the dual decomposition method 
determines that the combat units will not receive adequate 
logistical support, policies can be made to alleviate the 
problems of support. 

The stockage policy on the amount of supplies available 
at forward supply depots could be changed to meet the demands 
of the combat units. If the transportation network is 
inadequate, stockage levels could be increased to meet the 
demand. If the network is capable of carrying more supplies 
and assets are available, the stockage level at forward depots 
could be decreased to reduce the susceptibility of the depots 
to enemy forces. 

Policy concerning the amount of forces dedicated to the 
security of LOCs could also be made. The assignment of more 
security forces may decrease the probability of interdiction 
and thereby, increase the flow of supplies to the combat 
units. Security forces could also be distributed to different 
segments of the transportation network based on the analysis 


of the results determined by the dual decomposition method. 


B. FUTURE RESEARCH AREAS 
The dual decomposition method can be applied to other 


stochastic combinatorial optimization problems. These 
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problems find, from a finite set of alternatives, an 
alternative that optimizes the objective function. Examples 
of stochastic combinatorial optimization problems are 
sequencing, scheduling, and routing problems with random job 
durations or travel times. These problems require complete 
enumeration of all alternatives to find the optimal solution's 
distribution or moments. The dual decomposition method may 
be used to decrease the computational burden by limiting the 
number of alternatives that must be solved. This would allow 
combinatorial problems to be used on a more practical basis 
and without sophisticated computer resources. 

Variations on the four day logistics problem presented in 
this paper should also be explored. One variation is that if 
an arc has several different capacities that are dependent on 
the weather and the probability of each different weather 
factor is known, the dual decomposition method could be used 
to determine an approximation of the optimal objective 
function's distribution for a maximum flow problem. Other 
variations are that the demands could be random rather than 
known and the arc capacities could be from a continuous 
distribution rather than discrete. The dual decomposition 
method could be applied to these variations as well to provide 
an approximate distribution of the optimal objective function 


without full enumeration. 
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A major topic that requires further research is the 
determination of how a threshold value of the termination 
criteria should be established to provide a solution within 
a required tolerance level. This will enable the user of the 
dual decomposition method to be able to establish a threshold 


value that is consistent with the requirements. 


De) 


APPENDIX COMPUTER PROGRAM 


o sevivevesvke de lek kes Sesle ve feveve st ate e seve sesle tese es Yes ses ses ves fevesevesesves Seses esc sede Sese seve seve deve Se sekvveve sete slese te te Se Sesles Nese sie vey vetted sie 


THIS PROGRAM GENERATES THE FOUR DAY LOGISTICS PROBLEM ‘a NETWORK 
FORMAT AND USES THE DUAL DECOMPOSITION METHOD TO DETERMINE 
AN APPROXIMATE OPTIMAL OBJECTIVE FUNCTION 


* IT USES THE NETGEN, BUBSRT, NODPAK,AND HEDLST SUBROUTINES 
&: Sevevedevedevscesevetotoletetes sevevevedeve secede te sevetetecsetoleolose tee tet tose ke sete sek SISSON SONOS SESS ESSE: 


INTEGER. ‘MAXDEP MAXARC NXNARG Rees iey A noe, 


REAL DPINT( 50, O),PI PINTER 20, 50 atAR BC 190) Pret 300, 229300) 
INTEGER NN.NA.NDA.NIGM ND NX ARC NARC, TAIL(1 
CHEADC 100) ,CARL( 200, 1299)46 Pri 300 0 
CPITER( 10 42) EPNUM oe DEPNOD 700 sLEAP( 300), ,HCAP( 300), 
CDCAPL( 200,200) , DCAPH( 200,200) ,NH( 300); THRE 
CPAY( 3 Q) ANUM 300 
THTEGER NXT, MXDL,DUAL, LARC, TOT, PRICE,MXDH,HARC, PNTR 
*,..GNET INPUT DATA 
INTEGER*4 M, T( 300), wae C(300), CP(300), X(100) 
* PRIMAL, DUAL SOLUTIO 
INTEGER=4 Y(300), _UG190 
DARAMETER( MH NAR 00,MXNNOD = 100, MQLONG = 300,THRES = .5 ) 
INTEGER Dea MXNARC), -TAINNNOD = 100. Ur. TORL ror 
' CALI EXCNSC'F PILEDEF 08 DISK TH ANSWER Al (LRECL 130 PERM') 
MAXDEP = 100 
HAKAR Ce = 0 
' MAKDAY = 12 
CALL NETGEN 
* INPUT 
tot MAXDEP MAXARC ,MXNARC ,MAXDAY ,MXNNOD, 
9 NA,HN,NDA, 


NIGM, CAPL? CAPH, PINTER, DEPNUM, DCAPL,DCAPH,DPINT , DAY, ANUM, 
ARC; LCAP;HCAP,NH, PROB, 
* OUTPUT FOR’ GNET’ 

MNODRIEe) 


— NETWORK, CREATE ARC CAPACITIES 


NARC 
)G 0 1) THEN 


Ds, 
Bs bac 


CREATE RHS, INDUCE MAX FLOW FROM CAPACITIES TO LAST NODE 
DO L=1.MNOD 


ibe = 6 
Jz Aa): SED = 1 


= 1G) 


Ww 


END 
100 CONTIN 
NARC = 


xe 


Me x 
. 


peau) 


C 
D 
120 C 
x 


ne 


CON 


a SOT 


NARG = H(MNOD+1) 4 Or 


DSTI Me SS He Fe REPS SIS LTS he Se Seve seve deve se te sede dete seve ke ve sede deve de tek Sele Fe Fete Fe MO IT MM I Tete te 


54 


* SOLUTION OF THE MINIMUM COST NETWORK FLOW PROBLEM BY GNET 


2 

ate 

” 

tetokvedkdedededkededekdedes ve stetedevedeseds evede de deve dc devedesk ¥. ple ake ale ie tents 
Oe rreray 6S 30% 98 hy rie erly «6 35 Fe ray WS 9S FS VE ray WEDS OS FS BS UE OE 


BOTS ES BE tS Oh PR ae DE Te OB On ohn Din oo ee 


* THE DEST ARRAY Is THE QUEUE HOLDING THE PROBLEMS TO BE. “SOLVED 
MNOD | NARG)) 2 HjC.WGAr. xo y ®) 


CALL GNET( IpRT, IWwU 
DO H = 1-MNOD 
ii] u ri} 
o2 = RI + 0) - 1 
DO eeaLi(a) = UCIH) - UCT(J)) 
225 CONTINUE 
594 CONTINUE 
226) 1d =. 1. NARG 
PRICE = cd) “OY 3) 
TOTL = TOT + PRIC 
226 CONTINUE 
CAI GNET( IPRT,IWU, MNOD,NARC, T,H,C,HCAP, X, U,Y ) 
DO 324 IH = 1;>MNOD 
Jl = HCI) 
w2 = HiIH * 1) .- 1 
HO325 J =. 51,32 
DUALZ(3) = U(IH) - U(T(J)) 
325 CONTINUE 
324 CONTINUE 
DO 326 J = _1,NARC 
PRICE = Cty) % ¥¢J 
TO125= TO12 + PRIC 
326 CONTINUE 
Gun = 1 
Pe a al 
LF( i £705 Tor1 HON THRES } ) THEN 


TEC 08 oe oe TOT? TOT 


ELS 
CARI *PIVOT( DUALI, DUAL? . PNTR ,MENNOD.NARC Cc) 

CALL DESTOK( PNTR. DEST ,DUAL1. DUAL2 , NXT ,CUR,MXNNOD,TOT1,TOT2, 
INQLONG, PROB ,MXNARC,HCAP, CAB. NARC 


EWD 
ay CONTINUE 
CUR = N pA CUR .MQLONG , CUR) 


TFC CUR: EQ. NX 
ALL NOMORE 


Tuk Re eo 


DUALI(I BM ARG ey CUR I) 

LF( DEST Bul 21). ST ONT EN 
CP EST(CUR,I) 

Ss = HCAP(1) 


200 co UE 
Cue ee Ua ae MNOD,NARG, DT, HeG@ GR. x.) Way 1) 


“say 


HCE 
Se ee 

Bel os UCTH). = u(T(J)) 
J 


= ret), ‘ een 


ESTE CUR RC -2) 


ae Gea 
tl 


og 
Oo 
4 
Sh 
oun 
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DUAL2( I 


DEST( CUR 
Se I EN 


SCT! ONT 
2 dEstce Rea) 


= LCAP(I) 


R 


o~ I~ 


END I 
400 CONTINUE 
,»LWU, MNOD,NARG@, I,H,C,CP 9) Usr @) 


ee ZnoHnZs 
Q CRO CHOoSiwiH 
ce [el 


R a 
(tor ,TOT2,P,THRES) 
LT) HEN 


IRITE(4,*)TOT1,TOT2,P 


EN 
LF( | NOT. HLT THEN 
CALL PIVOT( DUAL] ,DUAL2, PNTR,MXNNOD.NARC) 
CALL DESTOR( PNTR. DEST, DUAL], DUAL2 NXT, CUR,MXNNOD, TOT1,TOT2, 
1MOLONG , PROB ,MXNARC ,HCAP, LCAP ,NARC) 


ELSE 

EUD IF 
GOTORTTE 
*. END safe ot 


., 0 9a ohn cba gti a¥a atain Ge ln cts abn cle alo abo cloatans arinla ate ale ale cloatr nda als aloaln als atniln cla clneke a0 clr ctnaUsinte Os cle at ecte ate ctocte lactate atnintnta iecty af Sate ataial 
she ste ste Sle le Ve le she she ste slo ale she a8e ale ah cle ale ale ae ake afe ale ale ale ale ale ale ale al ale ale ale ale ale ale ale ala ale alo ale ale ate Mo ale ale alo ale ale ale ale ale ole ale ale ale cle ale glo al ole ale ale ale ale ale ah 
DEBS BS GE GE BS BS BT BS OE GE BE BEBE BE GE OE BE SE DE IE FE BE BS BSE BE BS OS GE SE BE BE BE BE WE BE FES BE BS BS BE BE FE GE SE BS BE UE BE BE SE VE BE FE BE OE SE TS BE FE BE VE SE FE BE FE OE UE 


HLT_= 
IF 


go 
SOHM AsWESHaQ DoUWw 


Ing) 
" 
oO 
Oo > Mia 


* SUBROUTINE NETGEN 


ve ve ve Sete ve fe cle te deve ve ve te Seve ce ve ve ve sedis ve so ve ve ve dese seve te ve ve deste de ve Se eve seve dese ve ve ve See Se Se Fe Se Te se Le Fe Se Se Bee Pe ee See SE He Se 
J SUBROUTINE NETGEN 
we INPUT 

\ ( heir cetera SEIS Ao eee SES o NINO, 

3 NIGM, CAPL; CAPH, PINTER DEPNUM, DCAPL,DCAPH,DPINT,DAY,ANUN, 

4 NARC, LCAP ,HCAP,NH, PROB, 


OUTPUT. FOR’ GNET 
x MNOD,T,H,C ) 


aS aceite eae READS THE NETWORK AND CONSTRUCTS THE FOUR-DAY 


Me x4, 
; 


: vine ees te 

iy MAKARC = MAX ER OF ARCS 

7 MAXKDAY = MAX NUMBER OF DAYS 

ss MAXDEP = MAX NUMBER OF DEPOTS 

* MXNARC = MAX NUMBER OF ARCS IN TOTAL NETWORK 
s MXANNOD = MAX NUMBER OF NODES 

i OUTPUT ) 


; A = NUMBER OF ARCS 
* NN = NUMBER OF NODES 


* NDA = NUMBER OF DAYS 

* BIGM = A "LARGE" CAPACITY 

* CAPL = LOWER CAP OF ARC 

* CAPH = UPPER CAP OF ARC 
x PINTER = PROBABILITY OF INTERDICTION 

* DEPNUM = DEPOT NUMBER 

* DCAPL = LOWER CAP OF DEPOT 

* DCAPH = UPPER CAP OF DEPOT 

* DPINT = PROBABILITY OF INTERDICTION (DEPOT) 
* NH = HEAD 

ie SAIL 

* DAY = CORRESPONDING DAY FOR THE ARC 

* ANUM = ARCNUM FROM THE TABLE 
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* H = HEAD ENTRY-POINT ARRAY FOR GNET 


* INPUT 
& NUMBER OF NODES, NUMBER OF ARCS, NUMBER OF DAYS, BIGM 


* ARC #, LOWER CAP, HI CAP, PINTER 
a NUMBER OF DEPOTS 
us DEPOT NUMBER, NODE NUMBER, LOWER CAP, HIGHER CAP, DPINT 


REAL DPINT(MAXDEP ,MAXDAY ).PINTER(MAXARC, sMARDAY } ROB( 100) 
INTEGER NN,NA,.NDA.NIGM, NDPT.NXTARC,NARC, TAIL( 100), 
CHEAD( 100 ) peri 'beniopr| CAPH( 200 200i), 
CDEPNUM( XDEP) , DEPNOD 90) UMSNE. 
CDCA Lt 00.200 CAPH( 00, 00) NH(MXNARG 
CDAY(M NARC) oA AM MXN RC); LCA (85 ) HCAP( 300 

INTEGER LENGTH.C 400 al 0),A2(400),A3(400), 

A4( 400) , PRI PRIORC 3 0,300) ,SINK( 20) 


oo. GNELIOUIPUT D 
INTEGER*4 M, TA (MENARC), H(MXNNOD+1), C(MXNARC) 


CALL ExEMS FILEDEF 03 DISK NET DATA Al’) 
CALL EXCMS('FILEDEF 04 DISK TH OUTPUT Al (LRECL 130 PERM' ) 


READ(03,*)NN,NA,NDA,NIGM 
DO 1099 = 1, NA oA 
READ +) ARONUM TATL( ARCNUH) 4 READ REE oe) 

CEKPH ARGH NUM, J), PINTER( ARCNUM, J 
999 CONT 
1000 CONTINUE 
READ(O3, ,*)NDP 

ree a anger 


J 
03, S)DEPNUMCT), DEPNOD(1),DCAPL(1,J),DCAPH(1I,J),DPINT(I,J) 


998 RENT NUE 
1001 CONTINUE 
READ( 03 *) NUMSNK 
66 i £1, NN 
ice: = | FALSE. 
996 CONTINU 
DO 997 T = 1, NUMSNK 
SERED 2 4 
SNK 4 2 /TRUE 
IN ( ) =I 
ey CONTINU 
DO B80F daa t2NPA 
i+ (€J-1 = (10° TAIL(1)).+ J_- 
Ce Ret eeeh: A}) = E9 x TESS = @ = il 
NUCL * (C3; )).2 1 
bea +((J-1)*NA)) = J 
1004 CONTINUE 
1005 CONTINUE 


NXTARC (NDA*NA) + 1 


DO 1007 J = 
IMPNOD 5 INT TCI) LOD) 
1007 CONTINUE \ ie a mn oo 
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RC 
NETAEC =P 
Ax = oi SSEPRUM(Z) 
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NUMCNRTARG +t 
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ANUM NXTARG. © 
NXTARC = NXTARC + 2 
994 CONTINUE 
ELSE 


ENDIF 
995 CONTINUE 
NARC = NXTARC - 1 
CALL BUBSRT(MXNARC,NARC,NH, sNH DAY , ANUM) 
CALL NODPAK( MXNARC,NARC, ni ) 
DO 21 J=1,NARC 
1) CONTINUE 
CALL HEDLST 
1. ( MXNNOD,MNOD, MXNARC.NARG 
* OUTPUT ARRAY OF ENTRY POINTS BY HeAD’ NODE 


1 ee 
2 U 
Deo De ACR Soe re ere Oe ee See eae. Ae ay ae ey Daye os ecco re cas is as gms wae ae Coes rs eae Ca Ca as sas a 


e BUILD THE CAPACITY ARRAYS 


Fo sto we 5! fy ate x! fo alo at, ste by Po we lon! te 0 a! Wg cla Saic¥a aang ale abe cle clo ute ct foal, We stele ste whe Peale s’, 
Vevevedevedetevevetevedetetede vee ves es fe ve se le show! ate ate ale So ele 5 ve vede ves sleste seve seve veve ve se ve ve 
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DO 2000 page NARS 
IF( ¥CJ GT. ahi AND, ARUN 26% Q. AND. (DAY(J). LE. WA)) )THEN 
CAP GAL anon Ta) viGe 

HCAP = HCAN NUM( ds DAY( ) 

ELSE I ¢ AY} 0) THE 
LCAP = "NIGH 
HCAP(J) = NIGM 

ELSE IF(D ¥(J), ST, NA) THE 
LCAP( J) = CAPL(ANUM( J) ,DAY(NA 
HEAP(J) = CAPHCANUM(J ;DAY( NA 


END 
2000 CONTINUE 
bO 2 


O1 =1,NARC | 
IF( ANUNC de Wes ee L* ae 3) AY 1G 18, Oe 
LCA = BAe Gamo 
HCAP( J) = DCAPH( -1*ANUM DAY 
NEINU IF 
PEAR Rs Gohl ZENG ee oT epee aaa rises eels 


ate 
+ 


* READ IN DEMAND FOR EACH DAY 
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2209, CONTI 


* DO CAP ARRAY FOR DEMAND 
Sareircn ise 24 27 ace. ter 


DO 2001 NAR 
TE(SNK(-1 ~1 ANU UM 


2001 ote 
100 IT = i pau 


3 CONSTRUCT THE COST ARRAY 


sles slo ste oe ste te whe ate whe ate x e fe sie 5. 
nto oe Oe ey ve te ve deve we te ve de ve Seve se seve ve ve se seve ve ve 


‘DO 2210 I = 1,NDA 
2201 J = 1,NUMSNK 
READ(03,*)PRI 
PRIOR(I;J) = PRI 


E 

Ee 
a(t). 1 MA OvTHEN 

C = PRIOR(-1*DAY(1),SINK( ANUM(I))) 
G(1) = 0 


Sense cy cts oS. 
ste 


* DO PROB ARRAY 


Pa 


deve devs love docs Fede ve ve Ho deve te ves = 
DO 3000 
1F( (Day “shy on £0). AND. COR GT. 0), AND. (DAY(1). LE. NDA) ) THEN 


ELSE IF (baa), SE ee pA E) 2d. oy yrHEw 


ELSE TEC Rons(1},U0,0), Amp, (ANUM(T), GT. NDA) )THEN 


E 
PROB( 1) 
ND IF 
3000 Cae 
RETURN 
sersesesetey SUBROUTINE, BUBSRT(MANARC LENGTH CRIT AL Ag AB AS) Oo ccseaeseteeteaaaee 
“e 


* SUBROUTINE BUBSRT 


he 
ris 


dete de te de de vi te Xe ve de te de de de ve to de te te se de Be te te ve te te te de te ve te Be te te oe Fe Se Be te te Be Se te Me Be Fe ee Ne K Pete a STR TE ECE EE 
cANTHSRARE LENGTH, CRIT(MXNARC ) ,Al(MXNARC) ,A2(MXNARC) ,A3(MXNARC), 
- EGER See (00). DJ VCEME LD, TEMPZ. TEMPS, TEMPS, TEMPE 
DO sett J Sea eee 
COPLG@!) = CRIT(T) 
1000 CONTINUE 
DO OTGO1 1 = 1,LENCTH 
DO 1002 J_= I+1,LENGT 
DEC CORY( J )eLT, COPY( I) ) THEN 


59 


ie) 
© 
ing) 
tal 


ae 


Ina Il 


Ha 
Heo 
a 


Wo Te a he ta 


B 
Ns 
XY 
P 
P 
E 
P 
Al(I 
J 
i 
olf 
J 
J 
J 
J 


AAA > PP > PP} i | 
es]es] gales eel ey mre Ou) ey on 


EGE GEX 
WOO DGUUUGHHH 


> > > 
PWN PUN 
Woh 


ELSE 


JBINND)  1LJo 
1002 CONTINUE 
1001 CONTINUE 

RETURN 

Sener eee ie oe END se rede se vededesk seek Se See este ok ordeal ve ve Fe de Se See eA eee 
ate 
* SUBROUTINE HEDism 
severe? seve deve te ve ve de ve de ve ve ve ok de Sete ve oe Fe Ve te oe Bo Be Te oe Te Be Te Se ee Be Se BEE Ie Se EHS SOLE Se BEE TE TE SE TE TE OE 


SUBROUTINE HEDLST 
, 1 ( MXNNOD,MNOD. MXNARC,.NARC, 1T.NH, 
os OUTPUT ARRAY OF ENTRY POINTS BY HEAD’ NODE 


2 H 

y INTEGER MENNOD MNOD , MXNARC,NARC, T(MXNARC),NH(MXNARC), 
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